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Neutral stellar core at or over nuclear densities is described by a positive charged 
baryon core and negative charged electron fluid since they possess different masses 
and interactions. Based on a simplified model of a gravitationally collapsing or pul- 
sating baryon core, we approximately integrate the Einstein-Maxwell equations and 
the equations for the number and energy-momentum conservation of complete degen- 
erate electron fluid. We show possible electric processes that lead to the production 
of electron-positron pairs in the boundary of a baryon core and calculate the number 
and energy of electron-positron pairs. This can be relevant for understanding the 
energetic sources of supernovae and gamma-ray bursts. 

PACS numbers: 95.30.Sf, 87.19.1d , 23.20.Ra 

I. INTRODUCTION. 

In the gravitational collapse or pulsation of neutral stellar cores at densities comparable to 
the nuclear density, complex dynamical processes are expected to take place. These involve 
both macroscopic processes such as gravitational and hydrodynamical processes, as well as 
microscopic processes due to the strong and electroweak interactions. The time and length 
scales of macroscopic processes are much larger than those of the microscopic processes. 
Despite the existence of only a few exact solutions of Einstein's equations for simplified 
cases, macroscopic processes can be studied rather well by numerical algorithms. In both 
analytical solutions and numerical simulations it is rather difficult to simultaneously analyze 
both macroscopic and microscopic processes characterized by such different time and length 
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scales. In these approaches, microscopic processes are approximately treated as local and 
instantaneous processes that are effectively represented by a mo del- dependent parametrized 
equation of state (EOS). We call this approximate locality. 

Applying approximate locality to electric processes, as required by the charge conservation, 
one is led to local neutrality: positive and negative charge densities are exactly equal over all 
space and time. As a consequence, all electric fields and processes are eliminated. An internal 
electric field (charge separation) must be developed [l| in a totally neutral system of proton 
and electron fluids in a gravitational field. If the electric field (process) is weak (slow) enough, 
approximate locality is applicable. However, this should be seriously questioned when the 
electric field (process) is strong (rapid). For example, neutral stellar cores reach the nuclear 
density where positive charged baryons interact via the strong interaction while electrons 
do not, in addition to their widely different masses. As a result, their pressure, number, and 
energy density are described by different EOS, and a strong electric field (charge separation) 
on the baryon core surface is realized (2)] in an electrostatic equilibrium state. 

Furthermore, either gravitationally collapsing or pulsating of the baryon core leads to the 
dynamical evolution of electrons. As a consequence, the strong electric field dynamically 
evolves in space and time, and some electromagnetic processes can result if their reaction 
rates are rapid enough, for example, the electron-positron pair-production process of Sauter- 
Heisenberg-Euler-Schwinger (see the review [3)]) for electric fields E > E c = mlc 3 /(eh). If 
this indeed occurs, gravitational and pulsating energies of neutral stellar cores are converted 
into the observable energy of electron-positron pairs via the space and time evolution of 
electric fields. In this article, we present our studies of this possibility (the natural units 
h = c = 1 are adopted, unless otherwise specified). 

II. BASIC EQUATIONS FOR DYNAMICAL EVOLUTION. 

We attempt to study possible electric processes in the dynamical perturbations of neutral 
stellar cores. These dynamical perturbations can be caused by either the gravitational 
collapse or pulsation of neutral stellar cores. The basic equations are the Einstein-Maxwell 
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equations and those governing the particle number and energy-momentum conservation 

{n e , B Ul B ). v = 0, 

G^ u = —&ixG(T^ v + T^™), 

F^. v = AttJ", (1) 

in which the Einstein tensor G^ u , the electromagnetic field (satisfying F\ a p tl \ = 0) and 
its energy-momentum tensor T^™ appear; U" B and n e B are, respectively, the four velocities 
and proper number-densities of the electrons and baryons. The electric current density is 

J» = en p U^ - en e U», (2) 

where n p is the proper number-density of the positively charged baryons. The energy- 
momentum tensor T^ u = T£ u + is taken to be that of two simple perfect fluids repre- 
senting the electrons and the baryons, each of the form 

= p^gT + (Pe,B + Pe,B)U£ B U» B , (3) 

where p eiB (r,t) and p e>B (r,t) are the respective proper energy densities and pressures. 

In this article, baryons indicate hadrons, or their constituents (quarks) that carry baryon 
numbers. Electrons indicate all negatively charged leptons. Baryon fluid and electron fluid 
are separately described for the reason that in addition to baryons being much more massive 
than electrons, the EOS of baryons p B = Pb{pb) is very different from the electron one p e = 
peifie) due to the strong interaction. Therefore, in the dynamical perturbations of neutral 
stellar cores, one should not expect that the space-time evolution of number density, energy 
density, four velocity, and pressure of baryon fluid be identical to the space-time evolution 
of counterparts of electron fluid. The difference of space-time evolutions of two fluids results 
in the electric current ([2]) and field F^ v possibly leading to some electric processes. In a 
simplified model for the dynamical perturbations of neutral stellar cores, we approximately 
study possible electric processes by assuming that the equilibrium configurations of neutral 
stellar cores are initial configurations. 



i 



III. EQUILIBRIUM CONFIGURATIONS. 

In Refs. the equilibrium configurations of neutral stellar cores, whose densities are 
smaller than nuclear density n nuc i, are studied on the basis of hydrostatic dynamics of baryon 
and electron fluids in the presence of long-ranged gravitational and Coulomb forces. In these 
equilibrium configurations, very weak electric fields E <C E c are present, resulted from the 
balance between attractive gravitational force and repulsive Coulomb force. This electric 
field is too weak to make important electric processes, for example, electron-positron pair 
productions. We are interested in the case where strong electric fields are present. This 
leads us to consider strong electric fields in the surface layer of baryon cores of compact 
stars (quark or neutron stars) at or over the nuclear density. In this case, we assume that 
baryons form a rigid core of radius R c and density 
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where n p /fiB ~ N p /Nb < 1, Nb{N p ) is the number of total (charged) baryons and % 
?Vici ~ 1.4 x 10 38 cm~ 3 . The baryon core has a sharp boundary (r ~ R c ) of the width 
C ~ m~ l due to the strong interaction. The line element is Q] 

ds 2 = -g u dt 2 + g rr dr 2 + r 2 d9 2 + r 2 sin 2 9d(p 2 , (5) 
g -\r) = 1 - 2GM{r)/r + GQ 2 {r)/r 2 , 

where mass M(r), charge Q(r) and radial electric field E(r) = Q(r)/r 2 . 

Electrons form a complete degenerate fluid and their density n e e q (r) obeys the following 
Poisson equation and equilibrium condition j^]: 



d 2 V cq 
dr 2 



2 Id,, J dV r 



~ ~ n-rMdttgrr 



■ eq 



r 2 dr J dr 

= -AKeg rr (n p Ul-n c ?Ut), (6) 

E e = 9tt 2 \/\Pe\ 2 + rn 2 e - m e - eV cq = const. , 

where U p — XJ\ — (1, 0, 0, 0), , and Pf = (37r 2 ng q ) 1//3 are the Fermi energy and momentum, 
V eq (r) and E cq = —{g rr )~ x ^ 2 dV e(i {r)/dr are the static electric potential and field. In the 
ultrarelativistic case |P e F | 3> m e , we numerically integrate Eq. ([6]) with boundary conditions: 



ri e q ( r )|r«iJ c = n B 
i 

dr 
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As a result, we obtain on the baryon core boundary r m R c , the nontrivial charge-separation 
(n p — nl q )/riB and overcritical electric field E eq /E c > in a thin layer of a few electron 
Compton length A e [the curves (t = 0) in Fig. [Q. This is due to the sharpness boundary 
(C ~ m^ 1 ) °f the baryon core (HI) at t 
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Note that all electronic energy-levels 

^occupied = e I gl(, 2 drE eq (r) (8) 



are fully occupied and pair-production is not permitted due to Pauli blocking, although 
electric fields in the surface layer are over critical. We want to understand the space and 
time evolution of the electric field in this thin layer and its consequence in the dynamical 
perturbations of baryon cores, which can be caused by either the gravitational collapse or 
pulsation of baryon cores. 



IV. MODELING DYNAMICAL PERTURBATIONS OF BARYON CORES. 

It is rather difficult to solve the dynamical system (JTJ{3j) with the EOS ps = Pb{Pb) and 
Pe = Pe{pe) for the gravitational collapse or pulsation of baryon core and electron fluid, and 
to examine possible electromagnetic processes. The main difficulty comes from the fact that 
the time and length scales of gravitational and electromagnetic processes differ by many 
orders of magnitude. In order to gain some physical insight into the problem, we are bound 
to split the problem into three parts: (i) first, we adopt a simplified model to describe 
the dynamical perturbations of baryon cores; (ii) second, we examine how electron fluid 
responds to this dynamical perturbation of baryon cores; (iii) third, we check whether the 
resulted strong electric fields can lead to very rapid electromagnetic processes, for example, 
electron-positron pair production. 

As for the first part, we adopt the following simplified model. Suppose that at the 
time t — the baryon core is in the equilibrium configuration fl4j) with the radius R c 
and starts dynamical perturbations with an inward velocity R c (t) or pulsation frequency 
^puisa — R c /Rc- The rate of dynamical perturbations of baryon cores is defined as r~ u = 
R c /R c < c/R c . We further assume that in these dynamical perturbations, baryon cores are 
rigid, based on the argument that as the baryon core density n B)P (J4]) increases, the EOS of 
baryons Pb = Pb{Pb) due to the strong interaction is such that the baryon core profile (jl]) 
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and boundary width £ ~ m~ l are maintained in the nuclear relaxation rate T~ tTO ~ m n , which 
is much larger than r^. Thus, due to these properties of strong interaction, the dynamical 
perturbation of the baryon core induces an inward charged baryon current- density 

J B {R C ) = en p (R c )U r B (R c ), (9) 

on the sharp boundary of baryon core density (jl]) at R c , where the baryon density 
^B,p{Rc) = 0.5n Bp and the four-velocity U B (R C ) ^ 0. We have not yet been able, from the 
first principle of strong interaction theory, to derive this boundary property (Q of baryon 
cores undergoing dynamical perturbations, which essentially are assumptions in the present 
article, and the boundary density n BjP (R c ) and the boundary four- velocity U B (R C ) are two 
parameters depending on dynamical perturbations. This is in the same situation that so far 
one has not yet been able, from the first principle of strong interaction theory, to derive the 
sharp boundary profile (jl]) of baryon core densities of static compact stars [2[ . However, we 
have to point out that the boundary properties (jl]) and Q of the baryon core undergoing 
dynamical perturbations are rather technical assumptions for the following numerical calcu- 
lations of dynamical evolution of electron fluid and electric processes in the Compton time 
and length scales. These assumptions could be abandoned if we were able to simultaneously 
make numerical integration of differential equations for both dynamical perturbations of 
baryon cores at macroscopic length scale and strong and electric processes at microscopic 
length scale. 



V. DYNAMICAL EVOLUTION OF ELECTRON FLUID 



In this section, we attempt to examine how the electron fluid around the boundary layer 
of the baryon core responds to the dynamical perturbations of the baryon core described by 
the boundary properties (HI) and ([9]). Given these boundary properties at different values of 
baryon core radii R c , we describe electrons and electric fields around the boundary layer of 
baryon core by Maxwell's equations, the electron number and energy-momentum conserva- 
tion laws ([!]) in the external metric field (Q. In addition, we assume that the electron fluid 
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is completely degenerate, and its EOS is given by 



p e (t,r) = 2 I p°d 3 p/(2n) 
Jo 



where the single-particle spectrum is p° = (p 2 + m 2 ) 1 / 2 and the Fermi momentum is P e F = 
(37r 2 n e ) 1//3 . In the present article, for the sake of simplicity, we set the temperature of 
electron fluid to be zero and neglect all temperature effects, which may be important and 
will be studied in future. 

The electron fluid has four velocity Ug = (U t , U r ) e , radial velocity v e = (U r /U t ) e , U\ = 
Qtt^le and Lorentz factor j e = (1 + U r U r ) l J 2 = [1 + (g rr / ' gtt)^ 2 ]" 1 ^ 2 '■ 111 ^ ne res t frame at 
a given radius r it has the number density n e = n e y e , energy density e e = (p e + PeV 2 )le^ 
momentum density P e = (p e + p e )y 2 Ve 1 and v e = P e / (e e +p e ). In the rest frame, the number 
and energy-momentum conservation laws for the electron fluid, and Maxwell's equations are 
given by 

ne9tt 1/2 j + [neVe9tt 1/2 ) =0, (11) 



A 



(6e),t + (Pe), r + 
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dg rr dg tt ' 

——P e v e — (e e + Ve) 

dt dt y ' 



= -en e v e Eg tt 1/2 , (12) 

7-> 9rr \ , I - , n 9rr 

Pe + [Pe + PeV e 

9tt / t V 9tt / r 



. e e +p e ( dg tt dg rr 2 \ _ 1/2 

+ ^rUr-^J = - en ^ ' (13) 

(£),t = -4:ire(n p v p - n e Ve)g t t , ( 14 ) 

where (•••) ^ = ( — 9)~ 1 ^ 2 d(— 9) 1 {'")/9x, and in the line (|14p . the boundary velocity v p of the 
baryon core comes from the baryon current- density (J9J). We have the boundary four- velocity 
U B of the baryon core, 

u p = v fl = (CT/i7*)B, f/ B = ^ 1/2 7B, (15) 

and the Lorentz factor 

lB = (1 + £/ r £/^ 2 = [1 + (grr/M-W (16) 



s 



at the baryon core boundary R c . 

In the static case for v p = v e = 0, Eqs. fjl0til4|) are equivalent to Eq. (jSJ). Provided 
an initial equilibrium configuration ([6]) and proper boundary conditions, we numerically 
integrate these five equations ( 11011141) to obtain five variables n e (t,r),e e (t,r), P e (t,r),p(t,r) 
and E(t, r) describing the electric processes around the baryon core boundary. 



VI. OSCILLATIONS OF ELECTRON FLUID AND ELECTRIC FIELD. 

We consider the baryon core of mass M = 1OM and radius R c ~ 10 7 cm at the nuclear 
density n nuc i, and select its boundary velocity v p = 0.2c to represent possible dynamical 
perturbations of baryon cores. In the proper frame of a rest observer at the core radius R c , 
where gu(R c ) ~ 9rr(Rc), we chose the surface layer boundaries £_ ~ — A e ,£ + ~ 3.5A e , at 
which -Eeq(£±) ~ and proper thickness i = — and numerically integrate Eqs. ( [1011141) 
for the electron fluid. Numerical results are presented in Figs. [U and [2} showing that total 
electric field 

E(t,r) = E eq (r) + E(t,r), (17) 
where electron number density, energy density, and pressure oscillate around their equilib- 



rium configurations jj]. This is due to the fact that electrons do not possess the strong 
interaction and their mass is much smaller than the baryon one, as a result, the current den- 
sity of electron fluid in the boundary layer does not exactly follow the baryon core current 
density ((9]). Instead, triggered by the baryon core current ([9]), total electric fields E(t,r) 
deviate from E eq (r) and increase, which breaks the equilibrium condition ([6]), namely, the 
balance between pressure and electric force acting on electrons, dPf jdr + eE cq = 0. Accel- 
erated by increasing electric fields, electrons outside the core start to move inwards following 
the collapsing baryon core. This leads to the increase of the electron pressure (fTOj) and the 
decrease of the electric fields. On the contrary, increasing electron pressure pushes electrons 
backwards, and bounces them back. Overcritical electric fields work against the pressure of 
ultrarelativistic electrons. As a consequence, oscillations with frequency to = r~ ci ~ 1.5m e 
around the equilibrium configuration take place in a thin layer of a few Compton lengths 
around the boundary of baryon core. These are the main results presented in this article. 
We would like to point out that these results should not depend on the boundary properties 
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FIG. 1: The space and time evolution of the electric field (left) and charge-separation (right) 
around the boundary layer of the baryon core, M = lOM©, R c ~ 10 7 cm, and v p = 0.2c. The 
coordinate is £ = r — R c . 




FIG. 2: Time evolution of electric fields at different radial positions around the boundary layer of 
the baryon core, M = WMq, R c ~ 10 7 cm and v p = 0.2c. The coordinate is £ = r — R c . 



fllj) and (ED that we assume for the dynamical perturbations of baryon cores. The reason 
is that both electron and proton fluids in baryon cores are at or over nuclear density, and 
their Fermi momenta are the order of the pion mass m n ; therefore, electric fields must be 
at or over critical value E c = m 2 e /e to do work against motion of charge separation between 
positively charged baryon and electron fluids, and the frequency of oscillation because of the 
backreaction should also be the order of m e . It is worthwhile that these results are further 
checked by full numerical calculations without assuming the boundary properties fll]) and 
09]) of baryon cores, which undergo the dynamical perturbations caused by the gravitational 
collapse or pulsation. 

Suppose that the dynamical perturbation of the baryon core is caused by either the 
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gravitational collapse or pulsation of the baryon core, that gains the gravitational energy. 
Then, in this oscillating process, energy transforms from the dynamical perturbation of the 
baryon core to the electron fluid via an oscillating electric field. This can been seen from 
the energy conservation (TjQ) along a flow line of the electron fluid for v e ^ v p 

= en p F^U u B = en plelB (v p - v e )g rr E, (18) 

although we have not yet explicitly proved it. The energy density of the oscillating electric 
field is 

e osci = [E 2 (t,r)-E 2 cq (r)]/(8n). (19) 

The energy densities of the oscillating electric field and electron fluid are converted from 
one to another in the oscillating process with frequencies to ~ r~ ci ~ 1.5m e around the 
equilibrium configuration. However, the oscillating electron fluid has to relax to the new 
equilibrium configuration determined by Eqs. fl6]) and ((Tj) with a smaller baryon core radius 
R' c < R c . As a result, the oscillating electric field must damp out and its lifetime r re i ax is 
actually a relaxation time to the new equilibrium configuration. As shown in Fig. [2] the 
relaxation rate r~^ ax ~ 0.05m e . We notice very different time scales of strong interacting 
processes, electric interacting processes and dynamical perturbations of baryon cores: r^ ^> 

osci relax ^ coll " 

Moreover, when E(r,t) > E cq (r) (see Fig. [!)), the unoccupied electronic energy-level can 



be obtained by js] 



^unocuppicd 6 J 9rr drE(t,r^) ^ocuppied 

= e f gli 2 drE{t,r), (20) 



see Eq. (jSJ). This leads to pair production in strong electric fields and converts electric energy 
into the energy of electron-positron pairs, provided the pair-production rate r~^ T is faster 
than the oscillating frequency u = 7"~ ci . Otherwise, the energy of oscillating electric fields 
would completely be converted into the electrostatic Coulomb energy of the new equilibrium 
configuration of electron fluid, which cannot be not radiative. 
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VII. ELECTRON-POSITRON PAIR PRODUCTION 



We turn to the pair-production rate in spatially inhomogeneous and temporally os- 
cillating electric fields E(t,r). Although the oscillating frequency u is rather large, the 
pair-production rate rZ^ can be even larger due to the very strong electric fields E(t,r). 
The pair-production rate can be approximately calculated by the formula for static fields. 
The validity of this approximation is justified (see jj], 8]) by the adiabaticity parameter 
rj~ l = (u /m e )(E c / £ max ) <C 1, where E m£LX is the maximal value of the electric field on the 
baryon core surface r ~ R c . Therefore we adopt Eqs. (38) and (39) and (64)-(66) in Ref. 6| 
for the Sauter electric field to estimate the density of the pair-production rate in the proper 
frame at the core radius R c 

e 2 EE iEc/E)Go (a) e 2 EE 

^ 47T 3 G (<7) 47T3 ' {Zl) 

where E (instead of E) in the prefactor accounts for the unoccupied electric energy levels, 
G (a) — > and G (a) — > 1 for a — (£/X e )(E/E c ) 3> 1. The electron-positron pairs screen 
the oscillating field E so that the number of pairs can be estimated by A/" pa j r ~ ^R 2 C {E j e). 
The pair-production rate is r~ iT ~ 7^ pair (47ri?^)/A/'pair ~ am e (£/ \ e )(E / E c ) ~ 6.6m e > 
r~ ci . The number density of pairs is estimated by n pair ~ Af pair /(4:7iR 2 £). Assuming the 
energy density e OSC j of oscillating fields is totally converted into the pair energy density, 
we have the pair mean energy e pair = e OS ci/^ P air- Using the parameters v p ~ 0.2c, R c ~ 
10 7 cm, and M = 1OM , we obtain e osci w 4.3 x 10 28 ergs/cm 3 , n pair ^ 1.1 x 10 33 /cm 3 , and 
epair ~ 24.5MeV. These estimates are preliminary without considering the efficiency of pair- 
productions, possible suppression due to strong magnetic fields, and possible enhancement 
due to finite temperature effect. 



VIII. GRAVITATIONAL COLLAPSE AND DYADOSPHERE 

Up to now, we have not discussed how the dynamical perturbations of baryon cores can 
be caused by either the gravitational collapse or pulsation of baryon cores. Actually, we have 
not been able to completely integrate the dynamical equations discussed in Sec. [Til for the 
reasons discussed in Sees. [TJ and HVl Nevertheless, we attempt to use the results of electric 
field oscillation and pair production obtained in Sees. [V] |VI] and IVHI to gain some physical 
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insight into what and how electric processes could possibly occur in the gravitational collapse 
of baryon cores. For this purpose and in order to do some quantitative calculations, we first 
model the gravitational collapse of baryon cores by the following assumptions: 

1. the gravitationally collapsing process is made of the sequence of events (in time) oc- 
curring at different radii R c of the baryon core; 

2. at each event the baryon core maintains its density profile and sharp boundary as 
described by Eqs. fl4j) and (jlj). 

The first assumption is based on the arguments that (i) in the electric processes discussed 
in Sec. I VI} the charge- mass ratio Q/M of the baryon core can possibly be approaching to 
1, then the collapse process of the baryon core is slowing down and its kinetic energy is 
vanishing because the attractive gravitational energy gained is mostly converted into the 
repulsive Coulomb energy of the baryon core; (ii) then this Coulomb energy can be possibly 
converted into the radiative energy of electron-positron pairs as discussed in Sec. I VI II and the 
baryon core restarts acceleration by gaining gravitational energy. We have already discussed 
the second assumption in Sees. IIVI and I VI I Here we want to emphasize that (i) the sharp 
boundary properties fl4]) and (Q in the second assumption are technically used in order to 
numerically calculate the dynamics of electron fluid in the thin shell around the baryon 
boundary (Sees. |V] I VI I and IVIip ; (ii) in the gravitational collapse or pulsation of neutral 
stellar cores at or over nuclear density, these sharp boundary properties (jlj) and (J9]) should 
be abandoned in a more realistic model of simultaneously integrating dynamical equations 
of electron and baryon fluids over the entire stellar core at macroscopic scales. This turns 
out to be much more complicated and we will focus on this study in the future. 

On the basis of these assumptions, the boundary velocity v p (R c ) ( TT5|) and boundary radius 
R c [or boundary density nB, P {R c ) ©] of the baryon core at or over the nuclear density are no 
longer independent parameters, instead they should be related by the gravitational collapse 
equation of the baryon core. We adopt a simplified model for the gravitational co llap se of 
the baryon core by approximately using the collapsing equation for a thin shell jsl, lOj 



fiV fdR, 



dt 



GM 



(22) 



where at different radii R c of the baryon core, we define the charge-mass ratio 

£<2 = Q C V(G^ 2 M) < 1; Q c « = R 2 c E eq , (23) 
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FIG. 3: The estimate of the core collapsing velocity v p = R c = dR c /dt at different collapsing radii 
R c for the baryon core of mass M = 10 M Q . 



and 

n = l-(M/2R c )(l + e Q ) 

F = 1 - (2M/R C ) + (Q C VR C ) 2 . (24) 

The collapsing Eq. (j22p for the collapsing velocity R c is based on the condition that at each 
collapsing radius R c , the shell starts to collapse from rest. As a result, using these Eqs. (|2"2"1 
[241) we describe the sequence of events in the gravitationally collapsing process in terms of 
the collapsing velocities v p = R c = dR c /dt defined by (fT5|) and (fT6l) at different collapsing 
radii R c of the baryon core, as shown in Fig. [3j Thus, at each event the induced inward 
charged baryon current-density (J9l) is given by 

J r B = en p U r B « en p (R c Vt/F), (25) 

as a function of the collapsing radius R c . The strength of this charged baryon current 
density (1251) depends also on the ratio of the charged baryon number and total baryon number 



(N p /Nb), which varies in the gravitational collapsing process because of the (3 processes |ll |. 
In this article, the (5 processes are not considered and the charged baryon (proton) number 
N p is constant; we select two values N p /Nb ~ 1/38 or N p /Nb ~ 1/380 for the charged 
baryon current density Eq. (125)) . The collapsing process rate is = R c /R c ^ c/R c . If 
the dynamical perturbation of the baryon core is caused by the gravitational core pulsation, 
the pulsation frequency can be expressed as w pu i sa ~ R c / R c = r co n. 

In the sequence of the gravitationally collapsing process, at each event characterized by 
[R c , v p (R c )}, we first solve Eqs. ()6]) and ([7]) of the equilibrium configuration to obtain the 
number density (nf 1 ) and electric field (E eq ) as the initial configuration of the electron fluid 
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FIG. 4: The energy (left) and number (right) densities of electron-positron pairs at selected values 
of collapsing radii R c for M = 1OM and N P /N B ~ 1/38 (upper); 1/380 (lower). We select Rf 3 * ~ 
10 7 cm so that ng ~ ^nucl- 



and electric field. Then, with this initial configuration we numerically solve the dynamical 
equations (I10H14P to obtain the dynamical evolution of electron fluid and electric field within 
the thin shell (a few Compton lengths) around the baryon core boundary, described by 
Eqs. PI and (j25|) . As a result, based on the analysis presented in Sec. IVHI we calculate 
the energy and number densities of the electron-positron pairs produced at each event in 
the sequence of the gravitationaliy collapsing process. These results are plotted in Figs. EH 
Limited by numerical methods, we cannot do calculations for smaller radii. 

In addition, at each event in the sequence of the gravitationaliy collapsing process, using 
the Gauss law, Q = R^E, we calculate the charge-mass ratio Q/M averaged over oscillations 
of electric fields, Q/M < 1 as shown in Fig. The averaged charge-mass ratio Q/M is not 
very small, rather about 0.4 (see Fig. [5]), implying the possible validity of the first assumption 
we made that the gravitational collapsing process is approximately made of a sequence of 
events. In principle, at Q/M = 1 the gravitational collapsing process should stop, whereas 
the gravitational collapsing process is continuous for Q/M = without considering electric 
interactions. 

It is clear that the ratio N p /Nb becomes larger, the charged baryon current density 
(JSJ) or fl25|) becomes larger, and all effects of electrical processes we discussed in Sees. N] 
IVII and I VHI become larger. As shown in Figs. HI for the ratio N p /Nb ~ 1/38, the energy 
density of electron-positron pairs is about 10 31 ergs/cm 3 , and the number density of electron- 
positron pairs is about 10 35,6 /cm 3 . The mean energy of electron-positron pairs is e pa ; r = 
e sci/%>air ~ 10-50 MeV. While, for the ratio N p /Nb ~ 1/380, the energy density of electron- 
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FIG. 5: The charge-mass ratio Q/M averaged over oscillations of electric fields is plotted at different 
collapsing radii R c for the baryon core of mass M = 10 M@. 



positron pairs is about 10 30 ergs/cm 3 , the number density of electron-positron pairs is about 
10 346 /cm 3 , and the mean energy of electron-positron pairs e pa i r = e osc i/n pa i r ~ 10-50 MeV 
does not change very much. 

It this article, it is an assumption that the gravitationally collapsing process is repre- 
sented by the sequence of events: the baryon core starts to collapse from rest by gaining 
gravitational energy, the increasing Coulomb energy results in decreasing kinetic energy 
and slowing down the collapse process, the electric processes discussed in Sees. IVj [VI] and 
IVHI convert the Coulomb energy into the radiative energy of electron-positron pairs, and 
as a result the baryon core restarts to accelerate the collapse process by further gaining 
gravitational energy. This indicates that in the gravitationally collapsing process, the grav- 
itational energy must be partly converted into the radiative energy of electron-positron 
pairs. However, we have not been able so far to calculate all processes with very different 
time and length scales from one event to another in the sequence, so that it is impossible 
to quantitatively obtain the rate of the conversion of the gravitational energy to the energy 
of electron-positron pairs. Nevertheless, by summing over all events in the sequence of the 
gravitationally collapsing process, we approximately estimate the total number and energy of 
electron-positron pairs produced in the range R c ~ 5 x 10 5 — 10 7 cm: 10 56 -10 57 and 10 52 -10 53 
erg for the ratio N p /N B « 1/38; 10 55 -10 56 and 10 51 -10 52 erg for the ratio N P /N B « 1/380. 
These electron-positron pairs undergo the plasma oscillation in strong electric fields and 



annihilate to photons to form a neutral plasma of photons and electron-positron pairs 



12|. 



This is reminiscent of a sphere of electron-positron pairs and photons, called a Dyadosphere 
that is supposed to be dynamically created during gravitational collapse in Refs. fl^ . 
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IX. SUMMARY AND REMARKS. 

In the simplified model for the baryon cores of neutral compact stars, we show possible 
electric processes for the production of electron-positron pairs within the thin shell (a few 
Compton lengths) around the boundary of baryon cores that undergo gravitationally col- 
lapsing or pulsating processes, depending on the balance between attractive gravitational 



energy and repulsive electric and internal energies (see the numerical results in Ref. 14j). 
This indicates a possible mechanism that the gravitational energy is converted into the 
energy of electron-positron pairs in either baryon core collapse or pulsation. 

In theory, this is a well-defined problem based on the Einstein-Maxwell equations, 
particle-number and energy- momentum conservation (|TJ)- (J3j) , and equations of states, as 
well as the Sauter-Heisenberg-Euler-Schwinger mechanism. However, in practice, it is a 
rather complicated problem that one has to deal with various interacting processes with 
very different time and length scales. The approach we adopt in this article is the adiabatic 
approximation: the interacting processes with very small rates are considered to be adiabatic 



processes in comparison with the interacting processes with very large rates [15[]. Therefore, 
we try to split the problem of rapid microscopic processes from the problem of slow macro- 
scopic processes, and focus on studying rapid microscopic processes in the background of 
adiabatic (slowly varying) macroscopic processes. The adiabatic approximation we adopted 
here is self-consistently and quantitatively justified by process rates 

T strong ^ r pair > r osci ^ T rclax ^ T coll' (^6) 

studied in this article. In addition to the adiabatic approximation, we have not considered in 
this over simplified model the hydrodynamical evolution of baryon cores, the back-reaction of 
oscillations and pair-production on the collapsing or pulsating processes, and the dynamical 
evolution of the electron-positron pairs and photons. Needless to say, these results should be 
further checked by numerical algorithms integrating the full Einstein-Maxwell equations and 
proper EOS of particles in gravitational collapse. Nevertheless, the possible consequences of 
these electromagnetic processes discussed in this article are definitely interesting and could 
be possibly relevant and important for understanding energetic sources of supernovae and 
gamma-ray bursts. 
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